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Abstract. We discuss the nature of the two-stage percolation transition on the 
enhanced binary tree in order to explain the disagreement in the estimation of the 
second transition probability between the one in our recent paper ( J. Phys. A:Math. 
Theor. 42 (2009) 145001) and the one in the comment to it from Back, Minnhagen 
and Kim. We point out some reasons that the finite size scaling analysis used by them 
is not proper for the enhanced tree due to its nonamenablc nature, which is verified 
by some numerical results. 
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We have recently reported a numerical study of the two-stage bond percolation 
transition on the enhanced binary tree (EBT)p~|. Two percolation thresholds, p c \ ~ 
0.304 and p C 2 ~ 0.56, which respectively correspond to the divergence of the correlation 
mass and the correlation length, are obtained. The value of p C 2 estimated from the 
fractal exponent ip(jp) is consistent with the duality relation [2], p C 2 = 1 — p cl , where 
p cl ~ 0.436 is the lower threshold probability of the dual lattice of the EBT. On the 
other hand, Baek, Minnhagen and Kim estimated p& ~ 0.48 for the same model based 
on the finite size scaling (FSS) analysis [3]. This value is significantly smaller than 
our estimation while their estimation of p c i and p cl is consistent with ours. Thus they 
concluded that the duality relation does not hold for the EBT but inequality p C 2 < 1 — p cl 
is true. In this article, we compare these two estimations. In the following we use pb to 
note p c2 ~ 0.48 obtained in [3] for the distinction. 

First, we introduce the scenario of the second transition in the EBT, which has 
been already shown in pp. We only assume that connectedness function, Cq(£, p), which 
is the probability that a site at the £-th generation belongs to the same cluster with 
the root site, i.e., the site at 0-th generation, belongs to decays as a single exponential 
function; 

C (£,p) = A(p)2~^ = A{p)2^~ l)t , (1) 

at open bond probability p > p cl . Here £(p) is a correlation length and ip(p) = 1 — 

is a fractal exponent of the divergent clusters. We confirm the exponential decay of 

C (£,p) in Fig. [TJ Here we remarks on two quantities to detect the second transition, 

s Q {p,L)=Y,^Co{^p) and b(p, L) = 2 l ~ 1 Cq(L — l,p), (2) 

where L is a number of generations of finite size samples. We approximately identify 
x L — 1 with x L for x > 1 in the following, e.g., total number of nodes, N = 2 L — 1 — > 2 L . 
Substitution of eq. (CD) into eq. (TSJ) yields 

s °^ L ) = ^nrr iV " (p) and b ^ L ) = ^m N " iP) - ( 3 ) 

In these expressions, b(p,L) and so(p,L) are basically same quantities except 
unimportant coefficients and then we only treat b(p, L) in the following. Equation ([3]) 
leads to an important consequence that b is always infinite in the large size limit, 
N — ► oo, for p > p c i 0. Divergence of at p C 2, which is indicated in the right 
panel of Fig. [TJ results that ip(p) continuously approaches to unity to produce an 0(N) 
term [§|. What happens at p C 2 is essentially different from the ordinary second order 
transitions in amenable graphs. 

Next, we examine the analysis of Baek et al. in [3]. They assumed a FSS formula 

6(p,L)cxiV^3((p-^)iV 1/I7 ). (4) 
X The first threshold is defined by £(p c i) = 1 and then ip(p c i) = 0. 

§ Prefactor on Cq is possible but only results a correction factor (logTV) - ' 7 to sq and b. 
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Figure 1. (left) The connectedness function for six p's and four L's. Exponential 
decay can be observed before the boundary effect appears, (right) p-dependence of 
(the inverse of) the correlation length. Symbols indicate the values calculated by 
£{p,L) = — log 2 [Co(3L/4,p)/Co(i/4,p)]/[L/2] and dotted lines indicate the values 
calculated from 1 — ip(p) DP- The two estimation is almost same but the former 
is better near the p c \ to reproduce £(p c i) = 1- £ does not shows any singularity 
around p = 0.48 but approaches to zero at p w 0.56. (inset) The amplitude, 
A(p, L) = C a (L/2 lP )/2- L / 2 ^ L \ which hardly depends on p. 

This formula implies, in a sense of a standard FSS, that b is finite below pb and diverges 
as (pb —p)~^ v with infinite N. This seems strange because b has already diverged above 
p c i{< Pb)- Another diverging finite component which results a subleading term in b seems 
impossible since finite clusters growing with p must be absorbed to the already divergent 
clusters before diverges by themselves. We consider the scaling behavior is an artifact 
because eq. (J3J is approximately reproduced from eq. ([3]) without assuming another 
diverging component. Equation (J3J leads to b(p,L)/N* oc 2WW-Q L +*'(pb)(p-Pb)i+"-jf 
one chooses p& and <fi satisfying = ip(pb), b(p, L)/N^ looks a function of (p — Pb)L for 
\p ~ Pb\ *C 1 as 

b(p,L)^N%((p-p b )L). (5) 

This is obtained by replacing N^ 1 ^ with L = log 2 N in eq. (j4j). Note that L is locally 
approximated by a power function N 1 / u io^i(L) i7 local (L) = d In L/d In N = Lin 2, 
to reproduce eq. (jlj) in a narrow range of L. The two scalings are compared in Fig. [2J 
While the scaling with L shows good collapsing of data, the scaling with N l l u breaks 
down for large L (We use 1/V = 0.12 in [3] and treat larger generations by 7 than 
that in [3]) and only works in the narrow size range, around L = 12, as predicted from 
l/^iocai(12) ~ 0.120. Note that the scaling with L works for any pb G {pciiPdz) if 4> equals 
ip(p b ) (numerically confirmed too, not shown here) and therefore it does not gives the 
threshold of the second transition. Presumably some irrelevant finite size effect or short 
range behavior of Co yields the best FSS fitting point pb which depends on the data 
range of L. 

Another evidence for pb ~ 0.48 shown in [3] is the crossing of the ratio of the second 
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Figure 2. (left) Finite size scaling (FSS) corresponding to eq. (T4| using the parameters 
shown in 0; p b = 0.48, <p = 0.84 and 1/77 = 0.12. (right) FSS corresponding to eq. © 
using pf, = 0.48, <fi = 0.84. We show guide lines proportional to 2 -3 -°( p-Pi, ) L with light 
gray color. In both scalings, we use the Monte-Carlo data for 0.405 < p < 0.475 ( 
0.005 step ) averaged with 160000 samples. We show the same FSS of s& together. 



largest cluster to the largest cluster, (S2/S1}. Why crossing point gives critical point is 
based on the fact that the ratio (S2/S1) in the large size limit behaves as a step function 
of p around the critical point and takes a special value in the middle of step on the 
critical point, which is clearly confirmed by the FSS in the square lattice in [I]. Again 
it is not clear whether this is also true for the transition of the EBT. If the critical point 
between the non-percolating and percolating phases is replace by the critical phase, 
characterized by fractional ip(p), it is naturally expected that a slope appears to fill 
the gap. Such a slope is actually observed in the Cayley tree for p c \ < p < p C 2 = 1 
in [1]. Indeed we observe a tendency in the large L limit that (S2/S1) converges to a 
value which continuously decreases for p cl < p < p c2 rather than forms a step at pb (not 
shown here). In addition, we confirmed that (S2/S1) is far from a universal function of 
(p — PbjN 1 ^ (not shown here) unlike for the case of square lattice [I]. The crossing of 
($2/ s\) is considered to be caused by the change of the tendency in irrelevant finite size 
effect. 

In conclusion, we provided a simple scenario of the second percolation transition on 
the EBT and some numerical evidences which supports the scenario. We also showed 
that the FSS performed by Baek et al. does not holds for wide range of system sizes. 
Let us emphasize that the transitions of nonamenable graphs including the EBT is quite 
different from the usual second order transitions and standard analysis of second order 
transitions in amenable graphs cannot be applied directly to them. The value of p c2 
is, at least, larger than their estimation and the duality relation, p c2 = 1 — p c2 , seems 
valid for the percolation on the EBT. Baek et al. also claimed that the duality relation 
breaks down between the pair of {3,7} and {7,3} hyperbolic lattices based on the FSS 
analysis [3]. We also consider they underestimate the second threshold probability in 
this model. The duality relation should be true in this model since both of the dual 
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hyperbolic lattices are transitive in the large size limit [2] . 
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